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We present measurements of branching fractions in the b — > sss penguin-dominated decays 
B~^ — > 4>K~^ and B^ — > (j)K° in a sample of approximately 89 million BB pairs collected by 
the BaBar detector at the PEP-II asymmctric-cnergy _B-moson factory at SLAC. Wc determine 
B{B+ 4>K+) = (10.0^°;^ ±0.5) x 10"^ and ^ 0i^°) = {8At\l±0.5) x lO"''. Additionally, 

we measure the CP-violating charge asymmetry Acp{B^ — » (pK"^) = 0.04 ±0.09 ±0.01, with a 90% 
confidence-level interval of [—0.10, 0.18], and set an upper limit on the CKM- and color-suppressed 
decay B+ 07r+, B{B+ 07r+) < 0.41 x 10"® (at the 90% confidence level). 
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Decays of B mesons into charmless hadronic final 
states with a </> meson are dominated by & — s- sss glu- 
onic penguin diagrams (Fig. possibly with smaller 
contributions from electroweak penguin diagrams, while 
other Standard Model (SM) amplitudes are strongly 
suppressed Q. In the Standard Model, CP violation 
arises from a single complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix 2]. 
Since many scenarios of physics beyond the SM intro- 
duce additional diagrams with heavy particles in the pen- 
guin loops and new CP- violating phases a compar- 
ison of CP-violating observables with SM expectations 
is a sensitive probe for new physics. In the SM, ne- 
glecting CKM-suppressed contributions, the direct CP 
violation in P+ — > (f)K~^ 0|, detected as an asymme- 
try Acp = (r^K- - ^4,K+)/(r4,K- + ^4,K+) in the de- 
cay rates F^^^i = F(P^ (pK^), is expected to be 
zero; in the presence of large new-physics contributions 
to the b sss transition, it could be of order 1 [^. 
The B (j)K and P — > (/)7r decay rates are also sen- 
sitive to new physics; the latter is strongly suppressed 
in the SM, and a measurement of B{B — > (f)Tr) > 10~^ 
would serve as evidence for new physics Q . The branch- 
ing fractions of P"*" (pK^ and P° (pK^ have 
been studied by CLEO 0, PaPar HQ, and Belle [13 ; 
Acp{B+ 4)K+) has been studied by BaBAR 0. 




FIG. 1: Examples of quark-level diagrams for B — > (j>K and 
B (pn. Left: internal penguin diagram, right: flavor-singlet 
penguin diagram. 

This analysis is based on an integrated luminosity of 
about 82 fb~^, corresponding to approximately 89 mil- 
lion PP pairs, collected at SLAC with the PaPAR detec- 
tor at the PEP-II asymmetric-energy e+e~ storage 
ring operating on the T(4S') resonance. 

The asymmetric beam configuration provides a boost 
to the T{4S) in the laboratory frame {p-f w 0.56), in- 
creasing the maximum momentum of the P-meson decay 
products to 4.4 GeV/c. Charged particles are detected 
and their momenta measured by a combination of a sili- 
con vertex tracker (SVT), consisting of five double-sided 
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layers, and a 40-layer central drift chamber (DCH), both 
operating in a 1.5 T solenoidal magnetic field. The track- 
ing system covers 92% of the solid angle in the center- 
of-mass (CM) frame. The track-finding efficiency is, on 
average, (98± 1)% for momenta above 0.2 GeV/c and po- 
lar angles greater than 0.5 rad. Photons are detected by a 
CsI(Tl) electromagnetic calorimeter (EMC), which pro- 
vides excellent angular and energy resolution with high 
efficiency for energies above 20 MeV. 

Charged-particle identification is provided by measur- 
ing the average energy loss (dE/dx) in the two track- 
ing devices and by the novel internally reflecting ring- 
imaging Cherenkov detector (DIRC) covering the central 
region. A tt/K separation of better than 4cr is achieved 
for tracks with momenta below 3 GeV/c, decreasing to 
2.4(7 for the highest momenta arising from P"*" — > (f>h'^ 
decays. Electrons are identified with the use of the track- 
ing system and the EMC. 

We fully reconstruct P-meson candidates in the de- 
cay modes (j>h'^ and (pK*^, with <p K^K~ and K*^ 
TT+TT^ . For the track and the charged-track daughters 
of the (f> we require at least 12 measured DCH hits and 
a minimal transverse momentum px of 0.1 GeV/c. The 
tracks must originate from the interaction point (within 
10 cm along the beam direction and 1.5 cm in the trans- 
verse plane). Looser criteria are applied to tracks be- 
longing to Kg 7r~'"7r~. We combine pairs of oppo- 
sitely charged tracks originating from a common ver- 
tex to form K'^ and p candidates. A K'^ — > tt+tt^ 
candidate is accepted on the basis of requirements on 
the two-pion invariant mass (within 12 MeV/c^ of the 
nominal fsT^ mass 113), the flight-length (£) significance 
{£/(jg > 3), and the angle between the line connecting 
the P and K'^ decay vertices and the K'^ momentum 
(< 0.1 rad). Kaon tracks used to reconstruct the cp me- 
son are distinguished from pion and proton tracks using 
dE/dx information from the DCH in conjunction with 
dE/dx information from the SVT for track momenta be- 
low 0.7 GeV/c, and, for momenta above 0.7 GeV/c, with 
the measured Cherenkov angle and number of photons 
recorded by the DIRC. 

For an extended unbinned maximum-likelihood (ML) 
fit we parameterize the distributions of kinematic and 
topological variables for signal and background events in 
terms of probability density functions (PDFs). Each P 
candidate is characterized by the energy difference AE = 
{qr • qs/y/s) — and the beam-energy-substituted 

mass mES = [(s/2 + vr ■ Psf/E^ - |ll|. Here 

qr and qs are four-momenta of the T{AS) and the P 
candidate, s = (gr)^ is the square of the center-of-mass 
energy, pV and ps are the three-momenta of the T{AS) 
and the P in the laboratory frame, and Er = qj- is the 
energy of the T{AS) in the laboratory frame. For signal 
events, AP peaks at zero and mss peaks at the nominal 
P mass. The signal PDFs of both variables are ade- 
quately described by sums of two Gaussian distributions 
(whose means are not required to be the same). The 
background shape in AP is parametrized by a linear func- 
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tion and in m^js by a threshold function Candidates 
for our analysis are required to satisfy \AE\ < 0.2 GeV 
and mEs > 5.2 GeV/ c^. The variable AE provides addi- 
tional momentum-dependent tt/K separation in the ML 
fit for the — + branching fractions. The likelihood 
also incorporates the invariant mass of the 4> K^K^ 
candidate tiikk in the [0.99, 1.05] GeV/c^ range, which 
is described by a relativistic Breit-Wigner function con- 
volved with a Gaussian, a — l.OMeV/c^, determined in 
Monte Carlo (MC) simulation studies, to account for res- 
olution effects, and the cj) helicity angle 6*//, which is de- 
fined as the angle between the directions of the and 
the parent B in the rest frame. The cos Oh distribution 
is a quadratic function for pseudoscalar-vector B decay 
modes and is nearly uniform for the combinatorial back- 
ground. 

Backgrounds in the candidate sample arise primarily 
from random combinations of tracks produced in the 
quark-antiquark continuum. In such events, particles ap- 
pear bundled into jets, which can be identified with sev- 
eral variables computed in the CM frame. We use the an- 
gle 9t between the thrust axis of the B candidate and the 
thrust axis of the other charged and neutral particles ■ 
We require the angle 9t to satisfy | cos 6*^1 < 0.9. Other 
quantities that characterize the event topology are the 
CM angle 9b between the B momentum and the beam 
axis and the sum of the momenta pi of the other charged 
and neutral particles in the event weighted with Legen- 
dre polynomials Ln{0i),n = 0,2, where 9i is the angle 
between the momentum of particle i and the thrust axis 
of the B candidate. We combine these variables into a 
Fisher discriminant . Contamination from other B 

decays, as well as r+r~ and e~^e~jj production, is neg- 
ligible, as demonstrated in MC simulation studies. Pos- 
sible K^K~ S'-wave contributions, such as the /o(980) 
and the ao(980), are not expected to contribute under 
the (j) mass peak 1^ and are distinguished by their uni- 
form distribution in cos Oh', this systematic effect is small 
compared with current statistical and systematic uncer- 
tainties. 

We use an unbinned extended ML fit to extract signal 
yields and charge asymmetries simultaneously. The like- 
lihood for candidate j in the flavor category c is obtained 
by summing the product of event yield Nic and proba- 
bility Vic over signal and background hypotheses i. The 
total extended likelihood C for a sample of N events is 
given by 



N 



(1) 



The probabilities Vic are products of PDFs for each of the 
independent variables Xj = {mES, AE, T, rriKK, cos Oh}- 
The di are the parameters of the distributions in xj, 
which are fixed to values derived from signal MC, on- 
resonance sidebands in (toes, AE), and high-statistics 
data control channels i?+ 7r+Z?° (Z?° K^ir^) and 
B° TT+D- {D- K%-}. The control channels 



TABLE L Summary of branching fraction (B) and direct CP- 
asymmetry (Acp) results. A'sig and e are the signal yield and 
the total efficiency in the branching fraction fit. The 90% 
confidence-level interval for Acp is [—0.10,0.18]. 
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S (%) 


N ■ 


B (10-6) Acp 




6.7 




8A+\l±0.5 — 




19.6 


173 ± 15 


10.0lo;g±0.5 0.04 ±0.09 ±0.01 




20.4 


9+2-* 

"■^-0.9 


< 0.41 (90% CL) — 



have event topologies similar to those in cpK^ 
and i?'^ <t>K^, and are used to compare central values 
and resolutions of the variables toes, AE, and J- in data 
and MC simulation. By minimizing the quantity — ln£ 
in two separate fits, we determine the branching frac- 
tions, B, and the charge asymmetry, Acp, for (ph^ and 
4>K'^. In the (pKg case, there are two hypotheses, signal 
and background (i = 1,2), and a single flavor category. 
In the fit for B^ (j)h^ decays, we determine the flavor 
of the high- momentum track by comparing the measured 
Cherenkov angle with that expected for a pion or a kaon. 
In this way, the (ph^ {h = tt, K) decays are fitted si- 
multaneously with two signal {i — 1 for B^ — > (j)K^ and 
i = 2 for B^ — > (j)TT^) and two corresponding background 
{i — 3,4) hypotheses. We define the event yields riic in 
each of the two flavor categories (c = 1 for 5+ <j>h'^ 
and c ~ 2 for B^ —^ cph" ) in terms of the charge asymme- 
try Ai and the total event yield n^: nn = x (1 +Ai)l2 
and ni2 = x (1 — Ai)l2. 

For charged tracks originating from the interaction 
point, we determine the ratio of track-finding efficien- 
cies in data and MC simulation by conducting a study 
of a large sample of unambiguous charged-track candi- 
dates that have at least 10 measured hits in the SVT; 
the method relies on the fact that for both the SVT and 
the DCH the differences between the track-finding effi- 
ciencies in data and MC simulation are small, and so the 
two detectors can be used to calibrate each other. The ra- 
tio of TT+Tf- reconstruction efficiencies in data and 
MC simulation as a function of the K'^ momentum and 
decay point is determined from a study of a large inclu- 
sive sample of K'^ tt+tt" decays; this method employs 
the results of the tracking-efficiency study that covers K'^ 
decays occuring in the immediate vicinity of the interac- 
tion point. The charged- kaon-identification efficiencies 
in data and MC simulation are compared in a study of 
fully reconstructed D*^ D^'k^{D'^ K^ti^) decays. 

Results of the branching-fraction and CP-asymmetry 
fits are given in Table |l| Equal production rates of B^B^ 
and B^B^ are assumed. Figure|21shows the toes and Ai? 
distributions of (/)Xs(7r+7r-) and (f)K^ events together 
with the likelihood projections from the B fits. Goodness- 
of-fit tests have been performed to confirm that the values 
of likelihood C obtained in the fits are consistent with 
MC-based expectations. 

Systematic uncertainties in the ML fit originate from 
assumptions about the signal and background distribu- 
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FIG. 2: Projection plots of the variables rriES [(a) and (c)] and 
AE [(b) and (d)] in the fit for the (j)K+ (top) and (pKl ln+n-) 
(bottom) branching fractions. The data are shown by the 
histogram, while the curve is the result of the fit. The signal- 
to-background ratio is enhanced with a requirement on the 
signal probability 'Paig/(7-'sig -|- Pbkg) with the PDF for the 
variable being plotted excluded. 



tions and are dominated by the limited sideband and 
control-channel statistics. We simultaneously vary all 
PDF parameters within their uncertainties, and derive 
the associated systematic errors: 0.005 for Acp, 2.0% 
for B{(j)K+), and 2.8% for B{(j)K°). To account for the 
systematic uncertainty on the upper limit on S((/)7r+), we 
increase the upper limit by one standard deviation due to 
PDF variations (10.9%) and due to uncertainty in the re- 



construction efficiency (4.2%). The dominant systematic 
errors in the efficiency come from track finding (2.4% 
for B{(j)h+) and 4.2% for K(0iir°)), charged-kaon iden- 
tification (2% per 0), and reconstruction efficiency 
(2%). Other systematic errors from event-selection cri- 
teria, daughter branching fractions, MC statistics, BB 
backgrounds and _B-meson counting sum in quadrature to 
3.0%. The systematic uncertainty on Acp due to charge 
asymmetries in tracking and the DIRC is less than 0.01. 

In summary, we have studied branching fractions and 
charge asymmetries in the i?-meson final states and 
(i>K^\ the results are listed in Table U We do not observe 
a significant charge asymmetry in the mode B'^ — s- (pK^ 
and do not see evidence for B^ (pi:'^ . Our branching 
fraction and charge asymmetry measurements are consis- 
tent with, and supersede, our previous results reported 
in 13 . They are also consistent with existing SM pre- 
dictions. 
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